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Numerical Study on Aging Dynamics in the 3D Ising Spin-Glass Model. 
II. Quasi-Equilibrium Regime of Spin Auto-Correlation FunctionB 
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Using Monte Carlo simulations, we have studied isothermal aging of three-dimensional Ising 
spin-glass model focusing on quasi-equilibrium behavior of the spin auto-correlation function. 
Weak violation of the time translational invariance in the quasi-equilibrium regime is analyzed 
in terms of effective stiffness for droplet excitations in the presence of domain walls. Within the 
range of computational time window, we have confirmed that the effective stiffness follows the 
expected scaling behavior with respect to the characteristic length scales associated with droplet 
excitations and domain walls, whose growth law has been extracted from our simulated data. 
Implication of the results are discussed in relation to experimental works on ac susceptibilities. 
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§1. Introduction 

In recent years asiji£|di'ii£imics in spin glasses has been 
extensively studied.trBLlQB A basic experimental proto- 
col is isothermal aging which is relaxation after a spin- 
glass (SG) system is quenched from a high temperature 
above the SG temperature to a low temperature be- 
low Tc. Experimental studies of the relaxation of dc and 
ac susceptibilities during isothermal aging have revealed 
remarkable waiting time effects. The experiments have 
stimulated much theoretical interests and many theoret- 
ical ideas have been proposed most of which co 
count for qualitative features of the aging effects. 
Real tests on the theories are possible by investigating 
quantitative features of the aging effects, especially scal- 
ing properties. Progress in this direction, which may 
shed light on the controversial issues on nature of the 
SG phase, is still left io-he-done to a large extent. 

The droplet theoryja-Bli^l'llil' particularly the one pro- 
posed by Fisher and HusaS* contains concrete ansatz on 
possible scaling laws of the aging effects which are well 
amenable to be tested quantitatively by experiments and 
numerical simulations. In our previous paper (hereafter 
referred to as I) ,0 we have studied the relaxation of en- 
ergy of the three-dimensional (3D) Gaussian SG model 
with nearest neighbor interactions below using Monte 
Carlo (MC) simulations. We analyzed the data com- 
paring with the droplet theory. A key quantity in the 
droplet theory is the typical separation between domain 
walls R{t) which grows irreversibly by isothermal aging. 
Most of the scaling laws involve this characteristic length. 
Indeed in I we confirmed that the relaxation of excessive 
energy per spin Sexit) at time t with respect to the equi- 



* To appear in J. Phys. Soc. Jpn. 69 (2000) No. 4 
** Present address: Hydrographic Department, Maritime Safety 
Agency, 5-3-1 Tsukiji, Chuo-ku, Tokyo 104-0045 
*** E-mail: yhajime@ginnan.issp.u-tokyo.ac.jp 
•••• E-mail: takayama@issp.u-tokyo.ac.jp 



librium value follows the expected scaling law, 

6eT{t) ^ t(i?(t)/Lo)V(^WAo)', (1-1) 

when we used R(t) obtained by our independent simu- 
lation (see below). In the above equation T and Lq are 
characteristic energy and length scales respectively and 
d is the dimension of the space which is 3 here. The en- 
ergy exponent obtained is 9 = 0.20 ± 0.03 which agrees 
with the result of the defect energy analysis at T = 0.t3 
The simulated R{t) mentioned above was obtained 
by analyzing the replica-overlap functipBj— and found, 
in agreement with the previous worksE-illa) (but see 
alscO)), that the growth law of R{t) is well approximated 
by the form. 



(1.2) 



where tq is a certain characteristic unit of time. More 
explicitly, our data of R{t) at each T are well fitted to a 
function bt^^^ with adjustable parameters b and z. The 
latter depends on T and is expressed as 



l/z{T) ~ 0.17T, 



(1.3) 



at T < 0.7. Here and in the following, the Boltzmann 
constant fee is set to 1, and temperature T is measured 
by the unit of variance of the distribution of thfi interac- 
tion bonds J which is set to 1. As Kisker et alliS already 
pointed out, however, the simulated R{t) can be fitted 
also to a logarithmic function Rq + b{lnty/^ with ad- 
justable parameters i?o, b and ip, where ip is independent 
of temperature. This form, in the limit R{t) ^ Rq, is 
compatible with the asymptotic form proposed by the 
droplet theory. 



Rit) ^ Lo[iT/A)Ht/To)]'/^, 



(1.4) 



where A is a characteristic unit for energy barriers asso- 
ciated with thermal activation of droplet excitations. Lq 



and To in eqs.(1.2) and ( |l.4| ) are different sets of charac- 
teristic length and time scales. 
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The above-mentioned results may be summarized as 

1) Various quantities (such as 6eT{t)) simulated the 
isothermal aging and the characteristic length scales 
(such as R{t)) obtained by the independent sim- 
ulation satisfy the scaling forms (such as eq.(|l.l|)) 
which are predicted by the droplet theory. 

2) The growth law of R{t) is well fitted to the power- 
law form eq.(1.2), but it can be also fitted to a cer- 
tain logarithmic form which is compatible with the 



logarithmic growth law eq.(1.4). 
Feature 2) suggests that the computational time win- 
dow available by modern computers is not enough large 
for us to determine definitely a true asymptotic behav- 
ior of R{t) in the limit i — > oo. Therefore we call the 
above results seen in the present computational time 
window the pre- asymptotic behavior of aging. A simi- 
lar situation was already_encountered in the pioneering 
numerical work by Husell3^ who could confirm the den- 
sity of domain walls follows the expected scaling form 
1 / L'^{R{t) / LqY'^'^ with reasonable value of the fractal 
dimension ds but with the growth law being in some pre- 
asymptotic form. j— . 

In the present paper we study isothermal aging be- 
havior of the spin auto-correlation function in the so- 
called quasi-equilibrium regime. Here we focus our study 
on its scaling properties, or more explicitly, whether the 
simulated correlation function can be described by a cer- 
tain scaling form which is derived from the droplet the- 
ory, i.e., feature 1) above. To this end, we have per- 
formed extensive MC simulations on the 3D Gaussian 
Ising SG model in the same time window as that in our 
previous work I. 

The spin auto-correlation function is defined as. 



with 



Ciir; tw) = {Siir + t„)Si{t^)), 



(1.5) 



(1.6) 



where Si{t) is the sign of the Ising spin at site i at time 
t. At time t = 0, the system is prepared in a random 
initial configuration and we let it relax at T below Tc 
using conventional heat-bath MC dynamics afterwards. 
Here and in the following time t is measured by a MC 
step (MCS). By 1 MCS, all Ising spins are updated once. 
The auto-correlation function over some time intervals r 
at various waiting times <w are observed during a MC 
run of isothermal aging. The over-line in eq.(l.f;) de- 
notes the averages over sites and different realizations of 
interactions (samples), and the bracket in eq. (|l|) the 
average over thermal noises (or different MC runs) . 

The quasi-equilibrium regime, on which we focus in 
this study, is the time regime such that the time separa- 
tion T is relatively small compared with the waiting time 
Although behavior of the correlation function is not 
extremely different from the equilibrium limit within this 
regime, it exhibits systematic deviations from the ideal 
equilibrium behavior, namely, clear waiting time depen- 
dence, i.e., weak violation of time translational invari- 
ance. It is directly comparable with the measurement of 
relaxation of the Fourier spectrum of spontaneous mag- 



netic fluctuationstallS' in isothermal aging. Note that 
the Fourier component at a certain frequency lo is well 
defined only at large enough waiting times such that 
ojt^ 3> 1. For the same reason, the ac susceptibility 
is also considered to measure mainly relaxation of re- 
sponse of the system in this regime. From the point of 
view of the droplet theory, the quasi-equilibrium regime 
is defined more precisely as the time regime such that the 
typical size L{t) of droplet excitations which take place 
in the time scale of t is much smaller than the typical 
separation i?(iw) of domain walls present after waiting 
time On the other hand, let us recall for clarity that 
the so-called aging (or off-equilibrium) regime is the time 
regime such that r is as large as or much larger than 
the waiting time t^. More precisely, the aging regime 
is defined as the time regime where R{t = t + t„) is 
larger than R{tyj). The relaxation of the dc susceptibil- 
ities, which include the thermoremanent magnetization 
(TRM) and zero-field cooled susceptibilities (ZFC), can 
reflect strong non-stationarity within the aging regime. 

The droplet theory introduced a phenomenological 
concept called effective stiffness which characterizes re- 
duction of the excitation gap of small scale droplet ex- 
citations due to the presence of domain walls. The 
latter is considered to explain the waiting time effect 
within the quasi-equilibrium regime. We directly com- 
pare our results with this scaling ansatz. As far as we 
know, no direct quantitative analysis has been performed 
to test the scaling ansatz except for an experimental 
study of relaxation of Ihe ac susceptibility in a two- 
dimensional SG system. c£l' The previous numerical sim- 
ulations have analyzed the non-stationaiity_ p£ the spin 
auto-correlation within the aging regime,ErE2l'li3^ but the 
analysis of the weak non-stationarity within the quasi- 
equilibrium regime has not been carried out. 

We found that the spin auto-correlation function 
C{T;ty:,) in the quasi-equilibrium regime is consistently 
explained by the droplet theory in the following senses, 
i) As a function of L(r)/i?(iw), L{t) being the charac- 
teristic length of droplet excitations introduced above, 
it obeys a scaling law which is consistent with the sce- 
nario in terms of the effective stiffness-Q-* ii) Result i) is 
obtained when L{t) and R(t^) are assumed to obey the 
same growth law. The latter is a fundamental assump- 
tion within the droplet theory for spin glasses, while in 
a simple ferromagnet domain growth and droplet exci- 
tations are completely different kinds of processes, iii) 
We have confirmed result i) explicitly by making use of 
the growth law of R{t), which is obtained by our simu- 
lation and is approximated by the power law of eq.(1.2). 
Correspondingly, we have used 



L(r)^Lo(r/ro)i/^(^), 



(1.7) 



for L{t). These results add a further support for feature 
1) of the aging dynamics mentioned before, though they 
don't resolve the ambiguity pointed out as feature 2). 
We will argue that a similar situation to feature 2) is 
also encountered in the experimental studies up to now. 

The present paper is organized as follows. In the 
next section we briefly review the scaling ansatz of the 
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droplet theory. In §3 we present our simulated data on 
the spin correlation and response functions in isothermal 
aging. In particular, behavior of C(T;tw) in the quasi- 
equilibrium regime (r ^ t^) is examined in detail and 
is interpreted by the droplet theory. We discuss impli- 
cations of the present results with some other numerical 
and experimental results in §4. 

§2. Droplet Picture 

Here we review the droplet theory by Fisher-Huse 
(FH) concentrating on the influence of domain walls on 
the droplet excitations during isothermal aging. Let us 
begin by recalling the ideal equilibrium situation, where 
there should be no non-stationarity and th e be havior of 
the auto-correlation function defined in eq.(L5) becomes 
only a function of the time difference t. 

Equilibrium thermal fluctuations are considered as 
droplet excitationsliJ from a ground state, each of which 
is a global flip of a droplet (cluster) of spins within a 
distance L/2 from a certain given site i. The latter 
can be considered as a simple two-state system with a 
free-energy excitation gap ^l(?) and a thermal activa- 
tion time tl(z). The typical value F^f^ ^^"^ S^P ^^ii) 
scales as. 



T{L/Lof 



(2.1) 



where T is the stiffness constant of domain wall on the 
boundary of the droplet. The thermal activation time 
ri(z) is related to the energy barrier of the droplet BL{i) 
as, 

BLii)=Tln[TLii)/To]. (2.2) 

The typical energy barrier S^^^ of the energy barrier 
scales as, 

B'y^ ^AiL/Lof, (2.3) 

where exponent satisfies 9 < ip < d — 1. 

The probability distribution function pl{F) is 
sumed to follow a one-parameter scaling form. 



as- 



Pl{F) 




(2.4) 



where the typical gap f'^j^^ is given by eq.(2.1). A very 



important property of the scaling function p{x) is that 
it has finite intensity at a; = 0, p(0) > 0, which allows 
droplet excitations even at very low temperatures. The 
latter is considered to be responsible for ruany non-trivial 
properties of the SG phase in equilibrium^ It is also the 
case for dynamics in the quasi-equilibrium regime of the 
present interest. A similar one-parameter scaling form 
for ^'i(S) is assumed as. 



ptyp 



(2.5) 



where B^^^ is the typical energy barrier of eq.(2.3) and 



^'(a;) is a scaling function. 

The droplet excitations at different scales L which dif- 
fer by more than factor 2 are considered as approxi- 
mately independent from each other. Then using the 



two-state model, the equilibrium behavior of the spin 
auto-correlation function is obtained as, 



C(t) = lim C^{t-U 



Jl tanh^ 

L\TL{i)<T 



2T 



(2.6) 

where 11^ . . . represents multiplicative contributions of 
the relaxations of two-state systems at various scales L = 
2"Lo with n = 0, 1,2,...cxd enclosing the site i. The 
factor tanh^ (FL{i)/2T) is the equilibrium overlap of the 
two-state system of energy gap F^ii) with respect to the 
lower state (ground state). The constraint r^ii) < r is 
due to the fact that the droplets whose relaxation time 
is shorter than r can relax within the period of r while 
others with Ti(i) > t cannot relax appreciably and just 
give multiplicative contributions of 1. At low enough 
temperatures, the product in eq.( |2.6D can be expanded 
into a sum and the averages over sites (and equivalently 
over samples) can be performed using the probability 
distribution functions of the free-energy gaps and free- 
energy barriers. One obtains. 



1 - C{t) 

roo /.Tln(r/ro) 



La F Jo 



f°° dL f L 
where we introduced, 



dB^L{B) j dFpL{F)e-^'^ 
j dFpLiF)e-^/^, 
Pl{0), (2.7) 



^x) = / dy^{y), (2.8) 

and the characteristic length L{t) defined as, 

L(r)-io[(r/A)ln(r/To)]i/'^. (2.9) 

The function $(a;) satisfies $(0) = 1 due to the normal- 
ization of the distribution function 'i>{B) and is expected 
to remain close to 1 for x < 1 and decays rapidly for 
X ^ 1. The characteristic length L{t) is understood as 
a typical size of droplet which can flip in the period of r. 
In the last equation of eq.(2.7), we assumed that temper- 
ature is low enough so that only the term of order 0{T) 
contributes. 

Let us now consider the isothermal aging process. Dur- 
ing isothermal aging, there are domain walls separating 
different pure states, which are parallel or anti-parallel 
to a ground state, at a typical distance i?(tw) from each 
other. As far as one is monitoring short time ther- 
mal fluctuations due to small scale droplet excitations, 
such domain walls appear effectively frozen. The pres- 
ence of such frozen-in domain wall influences the small 
scale droplet excitations as the following. Some droplets 
which touch it can reduce their excitation gap compared 
with the others in the bulk of domains. FH considered 
the probability for such a circumstance to occur and the 
amount of reduction of the free energy gap of the droplet 
in such a circumstance. As a result, they found that in 



4 



Tatsuo KOMORI, Hajime YOSHINO and Hajime Takayama 



the presence of domain wall at a typical distance R from 
each other the free energy gap of a droplet of size L be- 
comes on average and typically as, 



Ffl^T''^[L/R]{LlLof 



L,R ~ ^ V^l ^^W-^l 1 (2-10) 
with the effective stiffness constant given by, 



T^ff[L/i?] = T ( 1 - c„ ( ^ 



(2.11) 



with Cv being a numerical constant. So as the system 
ages, it appears more and more stiff. 

Now we consider the behavior of the spin auto- 
correlation function in the presence of the domain walls. 
It can be obtained simply by replacing the typical free- 
energy excitation gap F*j^^ given in eq.(2.1) by the re- 



duced one F^j^^ given in eq.( 2.10| ). Particularly, we as- 



sume that the probability distribution of the free-energy 
gap FL ji[i) in the presence of the domain walls, which 
we denote as pl,r{F), follows the same scaling form as 
Pl{F) described in eq.(2.4). 



PL.AF) 



1 



ptyp 



ptyp 



(2.12) 



Replacing the distribution of the gap pl{F) in eq.(2.7) 
by Pl.r{F) defined above, we can obtain the modified be- 
havior of the spin auto-correlation function in the quasi- 
equilibrium regime L(t) <C R{t^) as. 



1 -C(T;iw) 



°° AL 



1 - Ceq(T) + 



L 

W) 

m 



T 



(2.13) 



T(L(r)/Lo) 
dy Hy) ( T 

Lo/L(r) y / 



1 - Ceq(r) + 



T-^{vL{t)IR{U)\ 
cip{0)T ( L{t) 



T{L{T)/Loy \Rit^) 



where the numerical constant Ci is given by. 



Cl 



(2.14) 



(2.15) 



In eq.(2.14), we expanded the effective stiffness given 
in eq.(2.11) assuming that L{T)/R{t^) is small enough 
in the quasi-equilibrium regime we are considering here. 
The last term is the major correction within the quasi- 
equilibrium regime to the equilibrium behavior described 
by the first two terms. The lower limit of the integra- 
tion Lo/L{t) is put to assuming L{t) is larger enough 
than Lq. Note that the integrand in eq.( p.l5 ) vanishes 
as ?/ — >• because of the inequality {d— l)/2 > 6x3 The 
integral is finite as far as the function <I>(.t) defined in 
eq.(^) decays fast enough for large x. 

The second term Coq(T) in eq.(2.14) iaJJie correlation 
function in equilibrium obtained by FH,E1P i.e., C(t) in 



eq.(2.7), 

Ceq(T) = lim C(r; tw) ~ <7ea + C2: 



Tp{0) 



■T(i(r)/io)^' 

(2.16) 

with qea being the EA order parameter given by 



<?EA = {Si)^ ~ 1 - Co 



(2.17) 



where q's are numerical constants. 

Let us next discuss briefly the relation between the 
correlation function C(t, t') with t — t + t„ and t' — 
above described and magnetic responses, or susceptibil- 
ities, frequently measured in experiments. In the quasi- 
equilibrium regime, even with weak violation of the time- 
translational invariance, the fluctuation-dissipation the- 
orem (FDT) in the following form is expected to hold 
at least approximately between C(i, t') and the response 
function G{t, t'), 



G{t,t') = 



1 dC{t,t') 



(2.18) 



T dt' 

From this equation the out-of-phase component of the ac 
susceptibility x"(w;tw) is evaluated as 

X"{uj;t^)^ / dtsin(wi) / dt'G{t,t')cos{ujt') 



2T 



(2.19) 



where C{Lo;tv,) is the Fourier component of C(T;t_). 
The latter is estimated, to a good approximation, agH3 



C(^;iw) 



TT 



-C{t;U) 



\lu\ dlnr 

with — 27r /lo. Hence we obtain 

TT d 



2T 91nr 



(2.20) 



(2.21) 



In the above derivation we have used the fact that 
C(T;tw), as a function of r and tw, varies quite slowly 
in a time scale of t^^. The fact that the FDT of 
eq.(2.19) holds well in the quasi-equilibrium regime, i.e., 
for R{t^) » L{t^), has already been tested by the care- 
ful experimental study of the spontaneous magnetic fluc- 
tuation which was directhi-compared with the relaxation 
of the ac susceptibility.llfl'lla-' The latter experiments have 
concluded that while there are obvious tw-dependences 
on both the magnetic fluctuations and ac susceptibility 
within the quasi-equilibrium regime, eq.(2.19) is satis- 
fied in the regime ut^ 3> 1 within the accuracy of the 
experiments. 

Substituting eqs.(??) and ( |2.16| ) into eq.( |2.2l] ), we 
obtain. 



x"(a.;tw)-XcqH ^ fHru) 



\Rit^) 



(2.22) 



where Xcq(^) is the equilibrium ac susceptibility given 
by 



m 



dL{T) I dlnrlT-^T^ 



(2.23) 



Numerical Study on Aging Dynamics in the 3D Ising Spin-Glass Model. II. 



5 



The above equation can be put into a simple form sug- 
gested by FHQ) using eg. (2. 11); 



with L{t) given by eq.(1.7) 



x"(w;<w) 



T<=ff[L(T^)/i?(t„)]-i - T- 



(2.24) 

This formula has been used in the experimental analy- 
sis of the. relaxation of the ac susceptibility in a 2D SG 
system.ES' 

Let us also note that from the FDT eq.( p.l8 ) we ob- 
tain the relation between C(i, t') and the zero-field cooled 
susceptibility XzFciT'jt^) as 



1 



dt'G{T + t^,t') = -[l-C(T;tw)]. 



XzFc(T;iw) = 

(2.25) 

Thus 1 — C(r;t„), which we will analyze in detail in 
the next section, is directly related to XzFc(''';iw) in the 
quasi-equilibrium regime. In the aging regime, on the 
other hand, the FDT no longer holds. The relation be- 
tween G(t, t') and C{t, t') in this regime has been exten- 
sively studied, mostly from the view-point of the mean- 
field theory,Ej^ but this problem is beyond a scope of the 
present work. 

To summarize, correction to the ideal-equilibrium 
behavior of the spin auto-correlation function is ob- 
tained based on the concept of the effective stiff- 
ness proposed by FH. The correction is a function of 
L(T)/R{t„) and remains small in the quasi-equilibrium 
regime L{T)/R{t„) ^ 1. Here L{t) is the typical size 
of droplet excitations which can flip within the time 
period r and R{t^) is the typical separation between 
the domain walls at waiting time t^. The time depen- 
dence of L{x) is given by eq.(2J). In the droplet the- 
ory by FH,S the growth law of i?(tw) is assumed to 
be the same (see eq.(1.4) ) based on the expectation 
that domain growth in spin glasses are driven by suc- 
cessive droplet excitations. On the time scale such that 
L{t) ~ R{tw), the above concept which separates the 
frozen-in domain walls and droplet excitations should 
break down, and the system crossovers to the aging 
regime ai R{t = t + t^;) ^ -R(tw) where much stronger 
non-stationarity is expected. 

Let us finish the review on the scaling theory with the 
following remarks. The above scaling form of C(r;fw) 
written in terms of i?(iw) and L{t) is expectecLtp hold 
more generally than the case described aboveO^ Par- 
ticularly, we expect it to hold for a certain class of the 
functional form of L{t) (and so R{tv,)), which is deter- 
mined by the probability distribution function ijB) in 
the above argument: a different l{B) f rom eq.(2.5) as- 
sociates with a different L{t) from eq.(2.9) which still 
ends up with the same scaling form of C(r;tw) written 
in terms of L{t) and R{t^). Actually, in our previous 
work I we found that ^' l (B) extracted from distribution 
of the largest relaxation time of small systems with sizes 
L has width which does not grow with L in contradiction 
to the scaling form of eq.( |2.5|) . However, this discrepancy 
is not crucial in deriving the scaling form of C(t; t^). In 
fact, from the restriction Ti(i) < t in eq.(2.6) alone, one 
finds the same expression for 1 — C{T;t^) as eq.(2.14) 



§3. Results of Simulations 

3.1 Model and method 

We have carried out MC simulations on isothermal 
aging phenomena in the same 3D Ising SG model as in 
our previous work I, i.e., the one with Gaussian nearest- 
neighbor interactions with zero mean and variance J = 1. 
The heat-bath MC method we use here is also the same 
as in I. The SG transition temperature is nunprically 
determined most recently as = 0.95 ± 0.04.ES The 
data we will discuss below are obtained dXT — 0.5 ^ 0.8 
in systems with linear size Lg — 24 averaged over 160 
samples with one MC run for each sample. In the pre- 
vious work I, it was confirmed that finite size effects do 
not appear for these parameters within our time window 
( < 2 X 10^ MCS). 

3.2 Spin auto- correlation function 

In Fig. 1, we show an overall profile of the raw data 
of the spin auto-correlation function C{T]t^) obtained 
by the present MC simulation. As can be seen in the 
figure, the curves of different waiting time t^ exhibit a 
characteristic crossover depeiuiing on t^ as already found 
in the previous studies.Bllil^cll' The data contain both 
the quasi-equilibrium regime and aging regime, and the 
crossover between the two regimes is expected to occur 
at around t ~ 
1 




10000 100000 



Fig. 1. The spin auto-correlation function C{T\t^) at T = 0.6 
as a function of t for fixed tw = 2(^"+'^) with n = 1 (tw = 
32), 2, 7{tw = 131072) from bottom to top, simulated on sys- 
tems with Ls = 24. 



According to the droplet theory, 1 — C{T:t^) in the 
quasi-equilibrium regime is given by eq.(2.14) with a con- 
tribution of the equilibrium part 1 — Ceq(r) and a cor- 
rection term which is proportional to {L{t) / R{t^)Y~^ . 
Using eq.(1.7) for L{t) in eq.(2.16), the equilibrium part 
is expected to behave as. 



Ccq(r)ocgEA + c(T/To)-^/^(^). 



(3.1) 



On the other hand, for thCjaging regime R{t = r + t^) » 
R{tvj) the droplet theoryQ^ predicts that C(r;tw) c>c 
{R{t^)/R{t = T + iw))^ where A is a p£w dynamical 
exponent which satisfies A > d/2. RiegerElP has pointed 
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out that with cq.( |l.2| ) for R{t), one obtains, 



(3.2) 



with A(T) — X/z{T). This interpretation has aUowed one 
to explain the apparent temperature dependence of the 
exponent A(r) of the data in the aging regiine_pbtained 
by MC simulations on 3D Ising SG models.EilllS For in- 
stance, using the values of exponent A(T) reported in 
Fig.2 of ref.E3 and our eq.(p^, A is evaluated at T < 0.7 
as 1.5 ^ 1.6 which is rather close to its lower bound 
d/2 = 3/2 = 1.5. 

Our interest in the present work is the details of 
C{T;tw) in the quasi-equilibrium regime. In order to 
investigate the latter regime, it is convenient to look 
at the spin auto-correlation function C(T;tv,) as a func- 
tion of iw with T considered as a parameter. In Fig. 2 
1 — C(T;t„) of some fixed small r's are plotted against 
tvj/r. Note that in this figure data at larger t^/r are 
closer to the equilibrium value as in the same sense in 
the experimental measurements of the ac susceptibilities 
during isothermal aging. According to the above crite- 
rion, the quasi-equilibrium regime is the part t^/r ^ 1 
where the droplet theory suggests scaling forms such as 
eqs.(2.14), (|^ and (|?2|). 
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and 



a[T) = 1 - Ceq(r), 



ci/5(0)| 



L{t)) \RiU) 



L{r) 



(3.4) 



(3.5) 



The correction term AC(r;tw) contains the factor 
{Lq/L(t)Y . For the data we have simulated this fac- 
tor is proportional to (t/tq)"^^^*^^-'. The te mpe rature 
dependence of exponent 1/ z{T) is given in eq.(l^). The 
energy exponent was obtained in our previous work I 
as 61 = 0.20 ± 0.03. Therefore the exponent e/z{T) be- 
comes quite small, and so the ratio of the above factor 
of the maximum time separation r and to that of the 
minimum r of the data shown in Fig. 2 is close to unity. 
Actually with Tmax/Tmin = 2^ and T = 0.6, the ratio is 
estimated as 0.94. This observation indicates that the 
factor {Lq/L{t)Y can be approximated as 



Lt(t) 



(7) 



e/z(T) 



1 



z{T) 



-In 



(3.6) 

where tq is a characteristic time scale corresponding to 
the length scale Lq, t,„ ~ (Tinin''max)^^^, and c a numeri- 
cal constant (~ (ro/rm)^/^^"^^). In eq.(^), in particular, 
it can be treated practically as a constant. 

Because of the above simplification, the correction 
term AC(r;tw) practically becomes a function of the 
scaled length L{T)/R{t^) = {t /t^f''^'^'^ alone. We now 
try to determine ^(t) for each r so that a master curve 
is obtained when AC(t;<:w) of different r are plotted 
against the rescaled time tw/^. Postponing the descrip- 
tion of details of the analysis till after eq.(3.8) below, 
we demonstrate a typical result in the inset of Fig. 2. As 
can be seen in the figure, this one-parameter scaling does 
work very well. 

0.25 I ' ■ 



0.2 



Fig. 2. The function 1 — C(T;iw) at T = 0.6 as a function of 
t„/T for fixed r = 2(2+") with n = 0, from bottom to 

top. The values of i^, for which we have simulated C(r;iw), 
are tw = 2^^°+"^^ with n = 1 (tw = 45), 2, 25{tw = 185360). 
When these data with each t are vertically shifted by the amount 
o(t) properly chosen, they lie top on each other as shown in the 
inset where AC(t; tw) = 1 — C{t; 4w) — <^ij)- 



0.1 



0.05 



10 



100 



1000 



3.3 Scaling analysis 

The scaling ansatz eq.(2.14) suggests that the data of 
1 — C(t; t^) can be decomposed into the equilibrium part, 
which does not depend on the waiting time and the 
correction term which depends on both t and t^. Let us 
denote the former as a(r) and the latter as AC(r;tw), 
i.e., 



1-C(r;iw) = AC(T;tw) + a(r), 



(3.3) 



Fig. 3. The extracted ^{t) of eq.(^) at T = 0.5, 0.6, 0.7 and 0.8 
from bottom to top. For T = 0.6 four sets of {air)} are drawn 
(see the text after eq. (3.8)). 



The chosen a(T) are drawn in Fig. 3, where we present 
also the result of the same analysis on C(r;tw) at T = 
0.5, 0.7 and 0.8 with l/z{T) = 0.085,0.115 and 0.127, re- 
spectively, as obtained in I. They appear as linear func- 
tions of Inr. This nearly logarithmic increase of a{T) can 
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be understood as the following, 
and (2.17) in §2, we obtain 



Combining eqs.(2.16) 



a(r) - 1 - Ceq(r) 

cTp{o) e 



z{T) 



T~p{0) 



ln(^ 



L{r) 



const. 



(3.7) 



where we used eq.(3.6) to derive the last expression. 
The above formula implies that slopes of a{T) with re- 
spect to ln(r) depend on temperature as because of 
the tempera ture dependence of the exponent l/z{T) as 
given in eq.(L3). From the data at four temperatures 
we obtained that da{T)/dhiT oc The stronger T- 

dependence than may be attributed to that of y5(0)/T 
in the last expression of eq.(3.7). 

It is noted that the second expression of eq.(3.7) im- 
plies that 1 — Coq(T) is concave in Inr as far as the EA 
order parameter ^ea is finite. But such a tendency can- 
not be seen in time-window of the present analysis. No 
tendency of saturation to gEA is seen also in Fig. 1 where 
an overall profile of the spin auto-correlation function ob- 
tained by the present simulation up to r < 2 x 10^ is 
shown. In fact up to now, no numerical data of the spin 
auto-correlation function which exhibit any tendency of 
saturation to some (/ea have been reported for the 3D 
SG raadeLs in spite of the enormous computational ef- 
forts.Eifcj) This is another pre-asymptotic behavior of 
the simulated data. 




L{t:)/R(Q={t/Q 



mT) 



Fig. 4. The double logarithmic plots of AC{T;t„) at T = 0.6 
versus L(T)/R{tm) with three different sets of {a{T)}: ai = 
o(t = 4) = 0.0690, 0.0702 and 0.0712 from top to bottom, which 
yield k ~ 2.2,3.0 and 3.3, respectively. 



Lastly let us discuss the exponent k defined by, 
AC(t;U) (L(r)/i?(tw))", 



(3.8) 



at the limit L{T)/R{t^) <^ 1. The questio n is if it 
agrees with the expected behavior given in eq.(3.5), i.e., 
K — d — 6. It is noted for clarity that the limita- 
tion of the computational time window does not al- 
low us to take equilibrium limit L(r)/i?(tw) — > ac- 
curately and thus leaves some ambiguities in the deter- 
mination of K. In Fig. 4, for examples, we show the 



double logarithmic plots of AC(T;iw) a,t T — 0.6 ver- 
sus L{T)/R{t^) — (r/iw)^^^^^^ for three different sets 
of {a{T)}. Here we used l/z{T = 0.6) ~ 0.102 using 
eq.(El). 



0.1 



O 0.01 




T=0.8 



r=o.6 



L{T)/R{Q={x/Q' 



Fig. 5. The double logarithmic plot of AC(t;4w) versus 
L{t)/ R{tY,). The slope of the dotted line is d (= 3). 



To draw these figures we have first chosen the value 
a{T = 4) (= ai) and then determined a(T) of the larger 
r subsequently in such a way that the data of AC(r; tw) 
lie on a single master curve. When we chose a too small 
value for ai, such as the case with ai = 0.0690 shown 
in the figure, we obtain a small a (~ 2.2) but a linear 
portion of the double logarithmic plot almost disappears. 
A too large ai (= 0.0712), on the other hand, ends up 
with a large n (~ 3.3) but the plot deviates from a linear 
behavior at smaller {L{t)/ R{tvj))- With a median value 
of ai (= 0.0702), which yields k ~ 3.0, we obtain linear 
behavior of the plot in a reasonable range of L{T)/R{t^). 
Thus the acceptable values for k from our simulated data 
are k 2.3 ~ 3.1. The latter values are consistent with 
the expected exponent in eq.(3.5), d — 9 with 9 ~ 0.20. 
However, it is hard to extract the value of by the data of 
AC(T;tw) alone. Let us remark however that the above 
ambiguity in k little affects the values of {a(T)}, in par- 
ticular, their slopes with respect to Inr are not sensitive. 
This is demonstrated in Fig. 3 by drawing four sets of 
{a{T)} at T = 0.6 with ai = 0.0690,0.0698,0.0702 and 
0.0712 which yield k ~ 2.2, 2.6, 3.0, and 3.3, respectively. 

These circumstances are similar at other tempera- 
tures we have examined. Therefore we show the rep- 
resentative results of {a(r)} in Fig. 3 which give rise 
to K ~ 2.8 ~ 3.0. In Fig. 5 we show AC(r;t„) at 
T = 0.6,0.8 with k ~ 3.0. For AC(T;iw) with a com- 
mon K, but with different temperatures Ti and T2, their 
relative magnitude in the range L{t) / R{tvj) ^ 1 can be 
estimated by that of {p{0)T/T)T=Ti through eq.(3.5). 
The relative magnitude of the two AC(r;tw) shown in 
Fig. 5 is about 1.3 which is roughly consistent with this 
estimation (= (0.8/0.6)1-4 ^ ^q^^ 

3.4 Out-of-phase linear ac susceptibility 

Let us finally discuss briefiy the logarithmic deriva- 
tive of the auto-correlation function, which is supposed 
to be almost identical to the out-of-phase linear ac sus- 
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ceptibility via FDT. Here let us denote the logarithmic 
derivative of the auto-correlation function simply as out- 
of-phase linear ac susceptibility x"(w;iw)- We show in 
Fig. 6 iw) at different lo which are obtained by nu- 

merical differentiation of eq.(2.21) using the same data 
of 1 — C(t; t^) shown in Fig. 2. 




Fig. 6. x"('^;<w) evaluated through eq.( |2.2l[ ) for = 2(2+") 
with n = 0, 1, ...,7 from bottom to top. In the inset the same 
data are plotted against tw/T^j. They lie on a universal curve. 



Using the growth laws eqs.(|1.2|) and ( pTT] ) in eqs.( p.22| ) 
and (2.23) respectively, the expected scaling behavior is 
obtained as, 



with 



oc 



{d-e)/z(T) 



(3.9) 



(3.10) 



where ujq = 2tt/tq. Note that the exponent 9/z{T) 
in eq.(3.1C ) originates from the factor {Lq/L{t)Y in 
eq.(2.23), and that the range of time r^^ = 27r/cLi of inter- 
est here is the same as that in the analysis of 1 — C(t; t^j). 
Because of the same reason as we noted above cq.(B.6), 
a master curve for x"(ti;;iw) has been obtained by sim- 
ply plotting them against /t^ as shown in the inset of 
Fig. 6. 

Let us remark that x"(^;iw) simulated here are ap- 
parently larger for the larger (smaller uj) . The scaling 
ansatz suggests that this is not the case for x"(w;tw) 
of much smaller frequencies, whose equilibrium values 
x"q(ci;) do depend on the frequencies and are the smaller 
for the smaller frequencies, in their asymptotic approach 
to the equilibrium values. We also note that, in order 
to check for the FDT to hold in the quasi-equilibrium 
regime it is desirable to perform direct simulation on the 
response against ac magnetic field. Since such an analy- 
sis is time consuming, it is left for a future work. 

§4. Discussions 

We have shown so far that the scaling properties of 
the simulated auto-correlation function C(r;<w) in the 
quasi-equilibrium range L{t) <C R{tvi) can be well ex- 
plained by the scaling ansatz by the droplet theory ex- 
pressed in terms of L{t) and R{t^). For the growth law 



of the latter, we have used the power-law form, eqs.( |1.2| ) 

More pre- 
d - 6*) in 



and (1.7) for i?(t„) and L{t), respectively, 
cisely, the fundamental exponents 6 and k (= 
the droplet theory, combined with the exponent \/z{T) 
of the growth-law of R{t^), consistently describe our 
simulated results: the energy relaxation in I and ag- 
ing behavior of C(t;<w) in the quasi-equilibrium regime 
in the present work. However, the time window of the 
present simulation and the previous one in I are in the 
range where the growth law of R.{t^) exhibits the pre- 
asymptotic behavior expressed as feature 2) in §1. Time 
scales of laboratory experiments on real spin glasses may 
also be in a similar range as we will discuss below. 

The quantity, which has been studied by experiments 
most frequently and is of interest from our present point 
of view, is the ac susceptibility x"(w;t). The latter is 
measured continuously in time t while the system is 
isothermally aged. Since the period = 27r/w is, in 
general, much smaller than the aging time t, it should be 
regarded as the response in the quasi-equilibrium range 
{t — + — iw) as we noted before. Indeed the scal- 



ing analysis based on eq.( 2.22 ) was directlK-jcarried out 
by Schins et al on a 2D Ising SG systcmEiJ^ Their re- 
sults axe consistent with the 2D version of the droplet 
theorycj) and, quite remarkably, the asymptotic growth 
law eq.(1.4) with a reasonable value of ^. 

Concerning the the ac susceptibility x"('^; in the 3D 
spin glasses, the Saclay group has often reported that 
x"(w;t) obey the cji-scaling; if x"{ijj\t) with different cj 
are vertically shifted by certain amount depending on 
UJ, and plotted against ut, they lie on a universal 
curveH'EZP The eguilibriui Ti valu es Xcq(^) thus obtained 
obey a power-law scaling, El H £3^ 



Xon(w) OC UJ . 



(4.1) 



The wt-scaling for ^("{uj^t) can be naturally explained 
by the scaling ansatz eq.( ^.22 ) with the power-law form 
for R{t) and L{t^), i.e., eq.(|3.9|). However, the data 
reported are not enough for examining the exponent 
k = k/z(T) = [d ~ e)/z{T) Th e above equation 
(4.1) just corresponds to our eq.(3.10) with a = 9/z{T). 
The time window of the experiments are < 10^^ and 
10^° < < 10^^ in unit of microscopic time Tc), 
while that of our simulation are ^ 10^ and 4 < < 
256. Presumably this difference allows one to extract the 
exponent a from the experimentally measured x"('^;i)- 
Intereaiijigly, a of CdIno.3Cr1.7S4 reported in Fig. 2b 
of Ref.Ea) is rather in good agreement with 9/z{T) eval- 
uated by our simulational studies on both the order of 
magnitude and the temperature dependence, though this 
is not the case at temperatures close to Tc nor at low 
temperatures T < 0.5Tc. In contrast to the above ex- 
perimental results which are favorable to the power-law 
form of R{tvj) and L{t), the Saclay group also reported 
that Xcqi^) observed in a wider range (7 decades) of uj 
fits rather well to .a-function of Inuj than to the power-law 



form of eq.( |3.10| ) .1 

Svedlindh et al^^ reported a Int-dependence of 
x"{uj',t) and a power-law dependence on cu in the log- 
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arithmic derivative of x"i^]t) as, 

dint ^ ' 



(4.2) 



at T < Tc (/3 ~ 0.25 at T = O.OSSTc)- This feature may 
also be explained by eg. ( 2.22 ) so long as L{t^) ob eys t he 
power-law form eq.(1.7). The exponent /? in eq.(4.2) is 
then evaluated as /? = {d — 26 )/ z, which is [{d/d) — 2] 
times larger than a in cq.( |3.10|) . 

From our point of view, another important quan- 
tity is the zero-field-cooled susceptibility XzFc{T;t„) or 
the thermoremanent magnetization t?t,trm(''"; tw) in the 
quasi-equilibrium regime. The two quantities are related 
by X7.FciT; t„) = [mpc - mTRM(T; t„)]/h, where mpc is 
the field-cooled magnetization under a field of strength 

and XzFciT'jtw) is directly related to 1 — C{T;t^) as 
expressed by eq.(2.25) in §2. Although xzFc(''';iw) (or 
'71trm(''"; ^w)) has been been measured fcequently since 
the pioneering work by Lundgren et al,c3 most of the 
measurements have been focused on its crossover be- 
havior between the quasi-equilibrium and aging regimes. 
Recently we have analyzed the data on mTRMiT',J^) of 
AgMn spin glass measured by the Saclay groupcZl^ by 
the same procedure as we did on 1 — C(t; t^) in §3. The 
results are qualitatively similar to those of x"{^',t) de- 
scribed just above: for the data at temperatures not close 
to Tc, the orders of magnitude of k and d as well as their 
temperature dependence are compatible with those ob- 
tained by our simulations, though this is not the case for 
the data at temperatures close to Tc- 

Besides the [experiments on the linear ac susceptibili- 
ties, Joh et altll' have recently extracted R{t) from the 
thermoremanent magnetization they have measured un- 
der fields of various strengths. Its growth-law agrees with 
the simulational result of eq.(1.2) even quantitatively, 
and the estimated R{t) with t of the order of laboratory 
time is at most only several tens of lattice distances, i.e., 
a length scale far from the thermodynamic limit we usu- 
ally imagine. 

Thus concerning the growth law of R{t) and L{t) in 
3D spin glasses, the experimental studies so far reported 
have not converg ed t o a unique conclusion: it is either a 
power-law as eq.(L2) or a logarithmic form as eq.(1.4). 
This ambiguous situation seems difficult to be resolved 
by looking at only u dependence of Xcq('^) because d in 
eq.(4T) is quite small partially because of the smallness 
of 9. The latter had already made it difficult to settle 
serious controversies concerning static properties of the 
SG phasca-* In this respect, it is of interest to exam- 
ine experimentally the cjt-dependence of x"(w;t) more 
in detail because /? in eq.( [4.2[ ) is relatively large. Also, 
as mentioned above, xzEciT", tw) in the quasi-equilibrium 
regime may be another important quantity to be studied 
in detail. Such experiments have to be performed in a 
time window of observations as wide as possible. 

Among problems on isothermal aging phenomena to be 
further explored numerically, a most important one may 
be aging dynamics at temperatures close to Tc, where a 
crossover from aging dynamics governed by the critical- 
ity at Tc to that governed by the T =^ criticality (the 
droplet theory) is expected to occur.lliP Indeed such a 



crossover in aging dynamics has been found xacently in 
the 4D Ising SG model by Hukushima et alE3> A sim- 
ilar analysis on the present 3D model may shed light 
on peculiar behavior of aging dynamics measured also 
experimentally at temperatures close to Tc as we have 
pointed out above. 

To conclude, we have simulated isothermal aging in 
the 3D Ising SG model, and analyzed behaviors of the 
spin auto-correlation function in the quasi-equilibrium 
regime. The simulated results are well interpreted by the 
scaling arguments which associate with the characteristic 
length scales i?(t„) and T(t), i.e, those of domains and 
droplet excitations in time scales of iw and r, respec- 
tively. We have conjectured that some aspects of the 
experimental data on the ac susceptibility and the zero- 
field-cooled susceptibility can be also explained along the 
same line. 
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